Single-File Diffusion in a Box 
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We study diffusion of (fluorescently) tagged hard-core interacting particles of finite size in a finite 
one-dimensional system. We find an exact analytical expression for the tagged particle probability 
density using a Bethe-ansatz, from which the mean square displacement is calculated. The analysis 
show the existence of three regimes of drastically different behavior for short, intermediate and large 
times. The results show excellent agreement with stochastic simulations (Gillespie algorithm). 



Introduction. - Recent advances in the manufacturing 
of nanofluidic devices allow studies of geometrically con- 
strained nano - sized particles in quasi one-dimensional 
systems in which crowding and exclusion effects are im- 
portant P, 0| ■ Situations where large molecules arc hin- 
dered to overtake also occurs in living systems as, for in- 
stance, protein diffusion along DNA [3[. Furthermore, bi- 
ological cells are characterized by a high degree of molec- 
ular crowding 

In this Letter, with experiments in mind, akin to 
we focus on diffusive motion of tagged finite-sized hard- 
core interacting particles (unable to overtake) (Fig. [1]). 
Such single-file systems show interesting behavior where 
the t x l 2 - scaling (t denotes time) of the mean square 
displacement S(t) = {[yr{t) — Vt,o] 2 ) in position yr{t) of 
the tagged particle [yr,o = 2/t(0)], for an infinite system 
with fixed concentration, is most striking ((.) denotes en- 
semble average). Also, the probability density function 
(PDF) p(yT,t\y T ,o) = pr is Gaussian @, 0, H- Even 
though single-file diffusion has received much attention 
H E3j El, to our knowledge very few exact re- 
sults are given for finite-sized particles in finite systems 
[l3j |. One exception is [l(| where the PDF for N dif- 
fusing point particles on a finite interval was obtained. 
However, asymptotic expressions were only given when 
the system was made infinite (keeping the concentration 
finite). Here, we go beyond previous results in the fol- 
lowing ways. First, finite-sized particles are considered 
and we show that the iV-particle PDF can be written 
as a Bethe-ansatz solution. Second, we perform a (non- 
standard) large N- analysis of pr, keeping the system size 
finite, showing, for the first time, the existence of three 
dynamical regimes: (i) t -c t co \\ = 1/ g 2 D where r co u de- 
notes mean collision time, D the diffusion constant, and 
g = N/ L particle concentration where L is the length of 
the system, (ii) r co u < i « r cq where r oq = L 2 /D is 
the equilibrium time, and (iii) t 3> r cq . Notably, only 
(i) and {ii) are found in infinite systems. Asymptotic ex- 
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FIG. 1: (Color online) Diffusing particles where mutual pas- 
sage is excluded, i.e yj < yj+i — A for j = 1, . . . , N — 1. 



pressions for pq- are derived in regimes (i)-(iii) and good 
agreement with (stochastic) Gillespie simulations [14]], is 
demonstrated. 

Statement of the problem. - We study N hard-core 
interacting particles with linear size A diffusing in a one- 
dimensional box of length L (Fig. [T]) . The probability of 
finding the particles at positions y = (y\, . . . , yjy) at time 
t, given that they initially were at yo = (yi,o, ■ ■ ■ ,yN,o), 
is contained in the iV-particle conditional PDF V(y,t\yn) 
which is governed by the diffusion equation 



dV(y,t\y ) 



D 



u 2 



i:) 2 



dt ' \dy\ ' ' ' dy 2 N 



V(y,t\y ). (1) 



Neighboring particles j and j + 1 are unable to overtake 
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dy 



3+1 



0, (2) 
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ensuring that yj+i — yj > A Vt provided that yjfi < 
2/7+1,0 — A. The boundaries at ±L/2 are reflecting 



D 



V(y,t\y ) 



dyi 



= D 

-(L-A)/2 



V(y,t\y ) 



dy 



N 



0, 



!/jv=(£-A)/2 



ambjorn@mit.edu 



(3) 

and the initial PDF is given by 

~P{y, °l2/o) = S(yi - 2/1,0) • • • S(y N - yN,o), (4) 

where 8(z) is the Dirac delta function. The tagged par- 
ticle PDF studied here is given by [l5| 

PT = — - — / dy[--- dy' N / dy' 10 ■ ■ ■ dy' N 
Pcq,ro Jr Jtzo 

x S(y T - y' r )8{y Tfi - y T . )V(y',t\yo)V e<l (yo) (5) 
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where p cq ,ra = j Uo dy' 10 ■ ■ ■ dy' N0 S{y T . Q - 2/^,0)^(^0 )> 
with integration regions 1Z — {yj+i — ])j > A,j = 
1, . . . , JV - 1; y x > —(L — A)/2; y N < (L - A)/2} and 
^0 = {Vj+1,0 ~ Vjfi > A,j = 1, . . . ,N - 1; 2/i.o > 
— (L - A)/2;y Nfi < (L - A)/2}. Initially, the particles 
are distributed according to the equilibrium density, 

= IJnj^fcft+i - w - A), (6) 

i.e the particles are distributed uniformly in the box. The 
function 9(z) is the Heaviside step function. 

Bethe - ansatz solution. - The (coordinate) Bethe- 
ansatz yields an integral representation of the PDF in 
momentum space 
Eqs. © - Q is 
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V(x,t\x ) 



The Bethe ansatz satisfying 

X • • • <j)(k N , X Nfi ) [ e ^ix l+ k^ + k 3 x 3 + ... + k N x N ) 

+ So je'^ 411 "'"' 11 '"'"'' 313 '""'"'"' 1 ''' 

where x±,..., xn and xi,o, ■ ■ ■ , xn,o are given in terms of 
y and yo according to (j = 1, . . . , N) 



= L - NA, x 



a ( ■ N + l 
Vj(fl) - A ( j — 



(8) 
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FIG. 2: (Color online) Particle trajectories, generated by the 
Gillespie algorithm, in a system where N = 11 . 



A — > 0, the iV-particlc PDF in [l(| is recovered. The 
single particle PDF is also found in [l(| (for A = 0) 
but in an unsuitable form for studies of finite systems. 
We obtained a more convenient expression, in which the 
large time limit is easily tractable, by finding the Laplace 
transform to ip(xj, xio; t), calculating the poles, and in- 
verting it back using the corresponding residues enclosed 
within the Bromwich contour 11511: 



where — 1/2 < x\ < X2 < . . . < xn < i/2 and 
— 1/2 < Xio < ... < xjv,o < 1/2. In fact, transforma- 
tion |(SJ) effectively maps Eqs. (HJ)- (|3J) onto a point-particle 
problem. The bracket in Eq. {7} contains N\ terms corre- 
sponding to all permutations of momenta k = k\, . . . , k^. 
The quantities Sij are scattering coefficients which con- 
tain information about the pair interaction between par- 
ticles I and j, and are in general functions of ki and 
kj . For the case of a pair interaction on the form given 
by Eq. ((2]), Sy = 1 (Sy = for non- interacting parti- 
cles) [l5j]. The time dependence enters through e - E i k ) 1 
with dispersion relation ("energy" ) E{k) = D (k\ + . . . + 
fcjy), obtained from Eq. (fTJ) . The functions (f)(kj,Xjfi) 
carry information about the boundary and initial con- 
ditions Eqs. ©-O, and for the finite box studied here 
<Kk j} x jl0 ) = 2^ l= - oc coa[k J (x^+i/2)}e^p+^, 
which was found using the method of images [15j. For 
an infinite system (£ — ► 00), cj)(kj,Xj t o) = e~ %k ' Xi -° [ill ]. 

Performing integrations over fci,...,fcjv, Eq. ||7J) is 
rewritten as 

V(x,t\x ) = ip(xi,Xi t0 ;t)ip(x2,X2,o;t) ■ ■ -il){x N ,x N $\t) 
+ ip(xi,X2 t0 ; t) ip(x2,xi )0 ; t) ■ ■ ■ ip(x N>0 , k N \ t) 
+ remaining permutations of £1,0, ■ ■ ■ , ^iv,0j (9) 



ie- Dk P is 



where tp(xj,xtfi;t) = J_ oo -^^(kj, x [fi )t 
the integral representation (inverse Fourier transform) of 
the (free) single particle PDF for particle j. Notably, as 



ip{xj,xi t0 ;t) = -|l + F m (t) 

x \v£> cos {— j^-J cos ^ — j — J 

. fxmxj\ . firmx lfi \-\ \ 
3m (^_^ sln (__^jjl )(1 0) 



-z/1 ^ sin I 



where = 1 ± (-l) m and F m (t) = e ^ {m7r)2J ^. 

Tagged particle density. - Integrating Eq. ((9|) accord- 
ing to Eq. ([5| leads to an exact form of the tagged 
particle PDF in terms of Jacobi polynomials P^'^\z) 
fl7j ] . given by [2(| 



N L \N R \ 
{{N l +Nr)iI>$(0,0,0;Q + N 



*(1,0,0;0 
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N R N L 



-*(0,1,0;0 

lpR.4> L 1pLlp R 



$(0,0, 



(ii) 



where £ = £(y T , yr.o, *) = ^r^l/^r^l), N l (N r ) de- 
notes the number of neighbors to the left (right) of the 
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tagged particle, and 

$M c-0 = (NL-(a + c))HN R -b)l NL - {a+c) 
9[a,t>>CU {NL + NR _ {a + b + c) yy V 



tc „{c,N R -N L +a-b) ( 1 + £ 
X ? r N L -(a+c) \l _ £ 



Arguments yx, UT,o an d t were left implicit. Also, 

-1 oo 



(12) 



^ = 
where 
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2 + T + ^ J m{x T ,x T fi)F m {t) (13) 
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Normalization gives ^ = 1 — t/>£ = ± — yj L 



(14) 



1 - ipi. and 



1 — -0j|, which completely determines p-r- A MATLAB 



implementation of pr using Eqs. pi[) - (fT4]) is available 
upon request. 

Figures [2] and [3] illustrate the typical behavior of the 
finite single-file system via stochastic simulations and 
pr- Figure [2] shows particle trajectories produced by the 
Gillespie algorithm (a Monte Carlo-like algorithm based 
on a lattice model which is equivalent to the master equa- 
tion (HI)- Figure [31 (a) illustrates the time evolution of 
pr for one tagged particle in the middle of the ensem- 
ble, and one by the edge. Snapshots of the PDFs are 
given at short (solid), intermediate (dashed) and large 
times (dotted). Notice the excellent agreement between 
the analytical result Eq. (fTTj) and the stochastic simu- 
lation. Panel (b) contains examples of the equilibrium 
PDF compared to the point-particle case A = 0. 

Three dynamical regimes. - Figure U depicts a nu- 
merical calculation of the mean square displacement 
of a tagged particle located in the middle of the en- 
semble. The solid blue lines were obtained from nu- 
merical integration (trapezoidal method) of S{t) = 



L/2-A(W H + l/2) 
L/2+A(Ar_ L + l/2) 



dyr (yr - yr,o) 2 pT, using Eq. (TT]) for 
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FIG. 3: (color online) (a) Tagged PDF [Eq. |TT}] depicted 
for the middle (red) and the rightmost (blue) particle, where 
NA/L = 0.25, at three instants of time. The initial condi- 
tions were y^fi/L = —0.15 and yiifi/L — 0.3. (b) Equilibrium 
density [Eq. (|16[l ] compared to the case of point-particles and 
a stochastic simulation (i/ r cq = 2) denoted by (o). The agree- 
ment between the simulations and the analytical results was 
checked using a x 2 -test with significance level a = 0.01. In the 
simulations, 500 lattice points and 10° ensembles were used. 



N = {3, 21, 141}. The behavior seen in Fig. g] illustrates 
the existence of three distinct regimes (i)-(m), which be- 
come more pronounced as N increases. 

In order to attain a deeper understanding of how 
regimes (i)-(iii) emerge, pr [Eq. pip ] was analyzed for 
large N, keeping L finite. A saddle-point approximation 
of the integral representation of $(a, b, c; £) [Eq. (fl~2j) ] 
proved unsuitable since it does not hold for all £ £ [0, 1] 
(i.e all times). However, making use of asymptotic forms 
of the Jacobi polynomial, derived in [l6l ]. we obtained a 
large iV-expansion of $(a, b, c; £) valid for all £ [2l[, and 
asymptotic expressions for pr in (i)-(m) and crossover 
times (t co ii and r eq ), were deduced: 

(i) Short times (t <C r co ii): For short times very 
few particle (wall) collisions have yet occurred and the 
particles are (to a good approximation) diffusing inde- 
pendently of each other. In this limit, pr is Gaussian 
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FIG. 4: (Color online) Mean square displacement for a tagged 
particle placed in the middle of the ensemble where N = 
{3, 21, 141}. Numerical calculations of S(t) based on Eq. (|11[) 
are indicated by thick blue lines. Straight lines show the ap- 
proximated analytic results for regimes (i) — {Hi), see text. 
Results from Gillespie simulations are denoted by (o) with 
errors 1 - S GaleBpie (t)/S(t) < 0.02. 



p T = (AnzDt)~ 1 ' 2 exp[-(y r - ?/T,o) 2 /( 40t )], with mean 
square displacement S(t) = 2Dt, which is in agreement 
with the numerical integration of Eq. (fTTj) . see Fig.0J 

(ii) Intermediate times (r co u <C t <C T eq ): The dynam- 
ics in this regime is dominated by particle collisions, and 
single- file behavior is observed: S(t) oc t 1 / 2 (Fig.0|). The 
PDF in this regime (for a particle located not too close 
to the edge) is found to be 



PT = 



2tt I 



4Dt 



1-gA 
B 



1/4 / 

exp 



(yr - vto) 



\ 



\ 



4Di 



1-gA 
Q 



the point-particle case [f| 0, Q to the finite particle case. 

(Hi) Large times (t ^> r oq ): For large times, pr reaches 
equilibrium [Fig. [3jb) contains examples], and S(t) is 
constant for t > r cq (Fig. |4|). The equilibrium density 
Pea.T is found using lim^oo £ = 1, leading to lim^i $ 
1 J22||, and a large t expansion of Eq. (| 13|) : 



(15) 



which is a Gaussian with a concentration dependent 
mean square displacement S(t) = [(1 — qA)/q] ^JADtpK. 
Thus, the simple rescaling g — > g/(l — gA) takes us from 



Peq,T 



1 



(N L + N R + 1)\ 



(L - NA) 



L 



N 



N L \N R \ 



x[-+y T - A(l/2 + N L ) 
[ | - y T - A(l/2 + N R ) 



Nr 



(16) 



Notably, Eq. (TIT))) is recovered by direct integration of 
Eq. ([6]), and also from simple entropy arguments [i~5l ]. 
The mean square displacement S(t — > oo 
case = iVfl, reads 



= iS eq for the 



5 nn — — 



JVh + 1 



L^VA\ 2 r(l/2 )r(2(jV fl + 1)) 
2 J T(N R 



i)r(JVfi 



5/2)' 
(17) 

where T(z) is the gamma function. 

Conclusions.- We have found an exact solution to a 
non-equilibrium many-body statistical mechanics prob- 
lem involving finite-sized particles diffusing in a finite 
system. The analysis showed, for the first time, the ex- 
istences of three distinctly different dynamical regimes 
for which exact analytical expressions of the PDF were 
found, using a non-standard asymptotic technique. The 
results showed excellent agreement with Gillespie simu- 
lations. 

The motion of tagged particles is sensitive to envi- 
ronmental conditions (e.g. concentration, diffusion con- 
stant and system size), suggesting that fluorescently 
tagged particles can function as probes or sensors at the 
nanoscale. 
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